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Abstract. The standard theory of metals, Fermi liquid theory, hinges on the key 
assumption that although the electrons interact, the low-energy excitation spectrum 
stands in a one-to-one correspondence with that of a non-interacting system. In the 
normal state of the copper-oxide high-temperature superconductors, drastic deviations 
from the Fermi liquid picture obtain, highlighted by a pseudogap, broad spectral 
features and T— linear resistivity. A successful theory in this context must confront 
the highly constraining scaling argument which establishes that all 4- Fermi interactions 
are irrelevant (except for pairing) at a Fermi surface. This argument lays plain that 
new low-energy degrees of freedom are necessary. This article focuses on the series 
of experiments on the copper-oxide superconductors which reveal that the number of 
low-energy addition states per electron per spin exceeds unity, in direct violation of the 
key Fermi liquid tenet. These experiments point to new degrees of freedom, not made 
out of the elemental excitations, as the key mechanism by which Fermi liquid theory 
breaks down in the cuprates. A recent theoretical advance which permits an explicit 
integration of the high energy scale in the standard model for the cuprates reveals the 
source of the new dynamical degrees of freedom at low energies, a charge 2e bosonic 
field which has nothing to do with pairing but rather represents the mixing with the 
high energy scales. We demonstrate explicitly that at half-filling, this new degree of 
freedom provides a dynamical mechanism for the generation of the charge gap and 
antiferromagnetism in the insulating phase. At finite doping, many of the anomalies of 
the normal state of the cuprates including the pseudogap, T— linear resistivity, and the 
mid-infrared band are reproduced. A possible route to superconductivity is explored. 



Mottness in High- Temperature Copper- Oxide Superconductors 



2 



1. Introduction 

Superconductivity in the copper-oxide ceramics stands as a grand challenge problem as 
its solution is fundamentally rooted in the physics of strong coupling. In such problems, 
traditional calculational schemes based on the properties of single free particles fail. 
Rather the physics of strong coupling resides in collective behaviour, signified typically 
by the emergence of new degrees of freedom at low energy. For example, in quantum- 
chromodynamics (QCD) the propagating degrees of freedom in the infrared (IR) are 
bound states not related straightforwardly to the ultra-violet (UV) scale physics. The 
key perspective presented here is that similar physics stems from the strong electron 
interactions in the copper-oxide superconductors. We will dehneate precisely how the 
emergence of collective behaviour at low energy accounts for many of the anomalous 
properties of the normal state of the cuprates. 

That the cuprates embody strong coupling physics stems from the Mott 
insulating[T] nature of the parent state. Such materials possess a half-filled band 
but insulate, nonetheless. Their insulating behaviour derives from the large on-site 
interaction two opposite-spin electrons encounter whenever they doubly occupy the same 
lattice site. For the cuprates[2], the on-site electron repulsion is typically U ~ ieV 
whereas the nearest-neighbour hopping matrix element is only t ~ 0.4eV. Although 
double occupancy is costly, there is no symmetry principle that forbids it even at half- 
filling. In the original proposal by Mott[T] to explain why NiO insulated, he assumed, 
as illustrated in Fig. ([T]), that each Ni atom remained neutral because 

U = E^+i + E^-i - 2E^ (1) 

dominates all other energy scales. Here is the ground state energy of an atom with 
valence electrons. For each Ni atom, N = 2. Hence, the zero-temperature state 
envisioned by Mott is one in which no atom is excited with an occupation of ± 1 
electrons. For NiO, this translates to no Ni+++ or Ni"*" ions exist as explicitly stated by 
Mottp. However, it is well-known[3l HI E] that the ground state of a Mott insulator 
possesses doubly occupied sites at half-filling. As a result, the simple cartoon [Ij that the 
Mott gap originates because double occupancy is forbidden is incomplete. Some have 
advocated [6] that in the Mott insulator, doubly occupied sites are immobile whereas in 
the metal they form a fluid. This account requires an explicit dynamical mechanism for 
the generation of the Mott gap. However, the dynamical degrees of freedom leading to 
the localization of double occupancy have not been unearthed. We offer here an explicit 
resolution of this problem. 

A few of the properties of doped Mott insulators are sketched in the phase diagram 
in Fig. ([2]). Aside from d^'^-y^ superconductivity, the pseudogap, in which the single- 
particle density of states is suppressed^ [H [13], and the strange metal, characterised 
by the ubiquitous T— linear resistivity [TH [15], stand out. As the phase diagram 
suggests, the pseudogap and strange metal phases are intimately related. That is, a 
correct theory of the pseudogap state of matter should at higher temperatures yield 
a metallic phase in which the resistivity scales as a linear function of temperature. 
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Figure 1. A half- filled band as envisioned by Mott. Each blue circle represents a 
neutral atom with N electrons and ground-state energy . The energy differences 
for electron removal and addition are explicitly shown. Mott reasoned that no doubly 
occupied sites exist because at zero temperature, U = + E^^^ — 2E^ ^ 0. This 

is, of course, not true. As a consequence the Mott gap must be thought of dynamically 
rather than statically. 

Nonetheless, numerous proposals fT7l \T8\ fT9| [20| EH [22] for the pseudogap abound 
that offer no resolution of T— linear resistivity. Part of the problem is that a series of 
associated phenomena, for example, incipient diamagnetism [23] indicative of incoherent 
pairing [201 [2ll [22l [21], electronic inhomogeneity[25l [261 [211 [SI [291 [IS], time-reversal 
symmetry breaking[30l[3ll|32ll33], and quantum oscillations [34] in the Hall conductivity, 
possibly associated with the emergence of closed electron (not hole) pockets in the first 
Brillouin zone (FEZ), obscure the efficient cause of the pseudogap and its continuity with 
the strange metal. Despite this range of phenomena, a key experimental measure [T5l [35] 
of the pseudogap onset is the temperature, T*, at which the first deviation from T— linear 
resistivity obtains. As a consequence, the physics underlying the strange metal must 
also yield a pseudogap at lower temperatures. Further, it must do so in a naturasl way. 
In our work, we take the relationship between the strange metal and pseudogap seriously 
and develop a theory[9l [101 E] that explains both simultaneously. In addition, we show 
that the same theory is capable of explaining other anomalies of the normal state such 
as 1) absence of quasiparticles[T2] in the normal state, 2) the mid-infrared band in the 
optical conductivity [361 [33 [3H1 [391 [101 [S], 3) spectral weight transfer across the Mott 
gap, and 4) the high and low-energy kinks in the electron removal spectrum. 

While it has been acknowledged for some time [12] that the normal state of the 
cuprates is incompatible with Fermi liquid theory, precisely what replaces it has not 
been settled. In a Fermi liquid, the low-energy excitation spectrum stands in a one-to- 
one correspondence with that of a non-interacting system. This correspondence must 
clearly break down in the normal state of the cuprates. The arguments of Polchinski [42] 
and others [131 HH [15] make it clear that breaking Fermi liquid theory in d = 2 requires 
new degrees of freedom at low energy, not simply 4-fermion interactions as they are all 
(except for pairing) irrelevant at the Fermi liquid fixed point. One possible origin of the 
new degrees of freedom[l6] is if spectral weight transfer between high and low energies 
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Figure 2. Heuristic phase diagram of the copper-oxide superconductors. In the 
strange metal, the resistivity is a linear function of temperature. In the pseudogap 
the single-particle density of states is suppressed without the onset of global phase 
coherence indicative of superconductivity. As discussed in Section 3.6, the dome-shape 
of the superconducting region with an optimal doping level of Xopt ~ 0.17 is not 
quantitatively accurate. See Fig. 19 for a more accurate determination of Xopt- 

mediates new electronic states at low energy. As a result, new states will emerge at low 
energy that have no counterpart in the non-interacting system. We show quite generally 
that this state of affairs obtains in the minimal model for a doped Mott insulator, namely 
the Hubbard model. Refinements of this model to include more details of the copper- 
oxide plane also retain this feature. We establish this result first through a simple 
physical argument which lays plain that in a doped Mott insulator, the phase space 
available for adding a particle exceeds the number of ways electrons can be added at 
low energy. Consequently, some new degrees of freedom not made out of the elemental 
excitations must reside in the low-energy spectrum. By explicitly integrating out the 
degrees of freedom far away from the chemical potential, the Wilsonian program for 
constructing a proper low-energy theory, we show that this new excitation is a charge 
2e bosonic field that in no way has anything to do with pairing. It is from this new 
degree of freedom that the pseudogap and T-linear resistivity follow immediately. Since 
this physics arises without any appeal to some further fact but relies only on the strong 
correlations of the doped Mott state, we have successfully isolated the efficient cause of 
the pseudogap. The associated phenomena mentioned above are supervenient on rather 
than central to the physics of the normal state. This review is organised as follows. In 
the next section, we discuss the experimental evidence for spectral weight transfer and 
show that it requires new degrees of freedom at low energy not made out of the elemental 
excitations. In section II, we derive the exact low-energy theory by formally integrating 
out the degrees of freedom far away from the chemical potential. In Section III, we 
compare the predictions of the theory with experiment. We close with a perspective on 
the remaining problem of superconductivity. 



Mottness in High- Temperature Copper- Oxide Superconductors 



5 



2. Mottness 

The origin of the Mott insulating state is subtle for two related reasons. First, the Mott 
gap cannot be easily deduced from the bare degrees of freedom in a model Hamiltonian. 
As remarked in the introduction, even in the Hubbard model, the ground state contains 
admixtures with the degrees of freedom, namely double occupancy, that lie above the 
gap. That is, if one were to write the bare electron operator[3] 

as a sum of two operators, one of which vanishes on doubly occupied sites, = 
(1 — ni_fj)cifj and its complement which is only non-zero when a site is doubly occupied, 
VitT = ni-fjCia, one would see immediately that such a separation is not canonical. 
As a result, T]i„ and C,ia have a non-zero overlap and hence they do not propagate 
independently as would be required for them to be gapped. In fact, it is unclear 
precisely how to write down a set of canonically defined fermionic operators that do 
become gapped as a result of the energy cost for double occupancy. This is the Mott 
problem. Its persistence has led LaughlinjlT] to assert that the Mott problem is fictitious 
and, in reality, does not exist. As mentioned in the preceding section, the real problem 
is that the Mott gap is fundamentally dynamical in nature. That it is difficult to write 
down the precise degrees of freedom that are becoming gapped is just a symptom of 
this fact. As will become clear from this review, the dynamical degrees of freedom that 
ultimately produce the Mott gap only appear when the high-energy scale is integrated 
out exactly. 

Second, all known Mott insulators order antiferromagnetically at sufficiently low 
temperature. To illustrate, two electrons on neighbouring sites with opposite spins can 
exchange their spins. This process proceeds through an intermediate state in which one 
of the sites is doubly occupied and hence the corresponding matrix element scales as 
t'^/U. Antiferromagnetism in the cuprates arises from this mechanism. This mechanism 
is distinct from the weak-coupling Slater [18] process in which a half-filled band orders 
as a result of nesting a.t Q = (7r,7r). While antiferromagnetism is certainly part of 
the Mott insulating story, it leaves much unexplained. It is that explanatory residue, 
namely the properties of Mott insulators which do not necessitate ordering, we refer to 
as Mottness. A simple property unexplained by ordering is the Mott gap itself. Above 
any temperature associated with ordering, an optical gap obtains p36t [371 Another 
such property is spectral weight transfer. 

2.1. Spectral Weight Transfer 

While what constitutes the minimal model for the cuprates can certainly be debated, 
it is clear [Hni EH [52] that regardless of the model, the largest energy scale arises from 
doubly occupying the copper dx2_y2 orbital. This orbital can hybridise with the in-plane 
Px and Py orbitals and hence a two-band model is natural. Since our emphasis is on the 
interplay between the high and low-energy scales, we simplify to a one-band Hubbard[3] 
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model 

^Hubb = - ^ I] 9ijclaCj,a + UJ2 ^T^i^^'i^^-T' (3) 
i,j,(T i,cr 

where i,j label lattice sites, gij is equal to one iff i,j are nearest neighbours, Ci^ 
annihilates an electron with spin a on lattice site i, t is the nearest-neighbour hopping 
matrix element and U the energy cost when two electrons doubly occupy the same site. 
Our conclusions carry over naturally to any n-band model of the cuprates as long as the 
largest energy scale is the on-site energy, U in Eq. ([3]) . That the dynamics of the charge 
carriers in the cuprates are captured by this model was confirmed by Oxygen Is x-ray 
absorption [53] on La2--a;Sra;Cu04. In such experiments, an electron is promoted from 
the core Is to an unoccupied level. The experimental observable is the fluorescence yield 
as a function of energy as electrons relax back to the valence states. The experiments. 
Fig- (E]), show that at a; = 0, all the available states lie at 530eV. As a function of 
doping, the intensity in the high-energy peak decreases and is transferred to states at 
528eV. In fact, the lower peak grows faster than 2x while the upper peak decreases faster 
than 1 — X. The separation between these two peaks is the optical gap in the parent 
insulating material. Though this observation of transfer of spectral weight from high to 
low energy is not expected in a semiconductor or a band insulator, it is certainly not an 
anomaly in strongly correlated systems. In fact, it is the fingerprint of Mottness as it 
has been observed in the classic Mott system NiO upon Li doping [19] and in all optical 
conductivity measurements on the cuprates [36l [371 [Ml [iQl HI] above any temperature 
having to do with ordering. 

This generic behaviour of spectral weight transfer is captured by the Hubbard 
model. To illustrate, consider the half-filled Hubbard model. A charge gap splits the 
spectrum into two parts, lower and upper Hubbard bands. Roughly, the lower Hubbard 
band LHB) describes particle motion on empty sites while particle motion on already 
singly occupied sites is captured by the upper Hubbard band (UHB). This relationship 
is only approximate because the UHB and LHB are mixed so that there are states in 
the LHB that have some doubly occupied character. To understand spectral weight 
transfer, we start in the atomic limit in which there is a clean gap of order U between 
the UHB and LHB. For a system containing N electrons on N sites, the weight of the 
LHB is N corresponding to N ways to remove an electron. The corresponding weight 
in the UHB is N as well as there are N ways to add an electron to the system. These 
bands are shown in Fig. Consequently, adding a hole in the atomic limit decreases 
the electron removal spectrum in the LHB by one state. The weight in the UHB is 
also affected as there are now N-1 ways to create a doubly occupied site. This leaves 
two states unaccounted for. The two extra states are part of the addition spectrum at 
low energies and correspond to the two ways of occupying the empty site by either a 
spin up or a spin down electron. In the atomic limit, the number of addition states 
scales as 2a; [Ml [55] when x holes are created. In a semiconductor or a Fermi liquid, the 
number of addition states would be strictly x. Experimentally [531 H9l [36| 1371 1391 HOlHT]. 
however, the low-energy spectral weight (LESW) grows faster than 2x. The excess of 2x 
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Figure 3. a) Normalized flourescencefSSj yield at the O K edge of La2-a;Sr2;Cu04+5. 
In the undoped sample, the only absorption occurs at 530eV, indicated by B. Upon 
doping the intensity at B is transferred to the feature at A, located at 528eV. b) 
Gaussian fits to the absorption features at A and B with the background subtracted. 
Reprinted from Chen, et al. Phys. Rev. Lett. 66, 104 (1991). 



can be understood simply by turning on the hopping[l]. When the hopping is non-zero, 
empty sites are created as a result of the creation of double occupancy. Such events 
increase the number of available states for particle addition and as a consequence the 
LESW increases faster than 2x. It is important to recall that the argument leading to 
the LESW exceeding 2x relies on the strong coupling limit. If this limit is not relevant 
to the ground state at a particular filling, the previous argument fails. 

2.2. Breakdown of Fermi Liquid Theory: More than just Electrons 

A natural question arises. Is spectral weight transfer important? A way of gauging 
importance is to determine if spectral weight transfer plays any role in a low-energy 
theory. A low-energy theory is properly considered to be natural if there are no relevant 
perturbations. Several years ago, Polchinski[42] and others [43l HH |45] considered Fermi 
liquid theory from the standpoint of renormalisation. They found [ISj HHl HH HH] that as 
long as one posits that the charge carriers are electrons, there are no relevant interactions 
(except for pairing) that destroy the Fermi liquid state. The setup [l2j is as follows. 
Decompose the momenta into the Fermi momentum and a component orthogonal to 
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Figure 4. Evolution of the single-particle density of states from half-filling to the 
one-hole limit in a doped Mott insulator described by the Hubbard model. Removal 
of an electron results in two empty states at low energy as opposed to one in the 
band- insulator limit. The key difference with the Fermi liquid is that the total weight 
spectral weight carried by the lower Hubbard band (analogue of the valence band in a 
Fermi liquid) is not a constant but a function of the filling. 



the Fermi surface 

p = k + l (4) 

Here 1 is the component orthogonal to the Fermi surface. Then consider scahng of 
energy and momentum towards the Fermi surface, in other words 

E^uE, k^k, e^ui, (5) 

where u is the scahng parameter. To quadratic order, the action is 

s = I dtd'p[tr{p)dMp) - {E{p) - EF{p))r ipmp)]. 

Hence, close to the Fermi surface 

E{p) - Ef(p) ^£vf, vf^ dpE (6) 
so that after scaling towards the Fermi surface (note that also t — > u~H) one finds that 

ip u-'^l'^i). (7) 

Consider now the four-fermion interaction. The argument to show that such interactions 
are irrelevant is particularly simple. In terms of powers of the scaling parameter, m, 
the measure over time contributes one negative power, the measure over the momenta 
orthogonal to the Fermi surface 4 powers and the 4-fermi interaction 4/2 negative 
powers. The delta function over the 4-momenta generically does not scale. Hence, the 
overall scaling of the four-Fermi interaction is governed by 1^-1+4-4/2 _ hence is 

irrelevant as the power of u is positive. The only exception to this argument if inversion 
symmetry is present is the Cooper pairing interaction. Consequently, as long as the 
charge carriers carry unit charge, there are no relevant interactions that destroy Fermi 
liquid theory. In the context of the cuprates, this argument is particularly powerful as 
it implies that in order to explain T— linear resistivity, some new emergent degrees of 
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freedom that have nothing to do with the electrons must be present. There have been 
attempts to circumvent this argument in the hterature that amount to essentially free 
field theory. In light of the above argument, such attempts must reduce to Fermi liquid 
theory and hence must yield resistivity. Others [56| [57] have directly confronted the 
Polchinski [42] argument and added extra derivative couplings to the Fermi liquid action. 
However, the relationship of such continuum models[56l [57] to any concrete realisation 
of Mott physics is not clear. 

It is straightforward to show that dynamical spectral weight transfer in a doped 
Mott insulator leads to a breakdown of the Fermi liquid picture and the emergence 
of new low-energy degrees of freedom. The interactions of the electrons with the new 
degrees of freedom can be formulated as a natural theory in which the electron spin-spin 
interaction is sub-dominant. As will be seen, the interactions with the new degrees of 
freedom govern all the physics that is independent of ordering. In this sense, we arrive 
at a natural separation between spin-ordering and Mott physics. To proceed, we define 
the number of single-particle addition states per site at low energy. 



as the integral of the single-particle density of states (A^(ti;)) from the chemical potential, 
/i, to a cutoff energy scale. A, demarcating the IR and UV scales. As long as A is chosen 
to exclude the high-energy scale, L is a well-defined quantity which simply counts the 
number of states in the unoccupied part of the spectrum at low energy. We compare 
this quantity to the number of ways an electron can be added to the holes created 
by the dopants. We call this quantity Uh- Our usage of 'ways' here refers to the spin 
degree of freedom of the electron only and not to combinatorics. From the perspective of 
single-particle physics, the intensity of a band is always equal to the number of electrons 
the band can hold. Hence, strict adherence to the single-particle picture requires that 
L = Uh, implying that the number of low-energy addition states per electron per spin 
is identically unity. For example, as shown in Fig. ([5]), in a non-interacting system, 
L = 2 — n = Uh- The same is true for a Fermi liquid as can be seen from the fact that 



Since the integral over all energies must yield 2, it follows that L = 2 — n. Hence, strictly 
for a Fermi liquid, L/ri/i = 1 as is dictated by the basic Landau tenet that the number 
of bare electrons at a given chemical potential equals the number of Fermi excitations 
(quasiparticles) in the interacting system. 

A doped Mott system is quite different because the total spectral weight in the 
lower-Hubbard band is not simply 2 but rather determined by the electron filling. 
Consider a Mott system in the atomic limit. As shown in the previous section, L = 2x 
is the exact result in the atomic limit because creating a hole leaves behind an empty 
site which can be occupied by either a spin-up or a spin-down electron. Likewise, the 
number of ways electrons can occupy the empty sites is rih = 2x. Hence, even in the 
atomic limit of a doped Mott insulator, L/nh = 1. However, real Mott systems are not 




(8) 



/ N{u)duj = n. 
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Figure 5. a) Integrated low-energy spectral weight, L, defined in Eq. ([5]), as a function 
of the electron filling, n: 1) the dashed line is the non- interacting limit, vanishing on- 
site interaction (f/ = 0), in which L — 2 — n, 2) atomic limit (blue line) of a doped Mott 
insulator, U = oo, in which L = 2(1 — n) = 2x, x the doping level and 3) a real Mott 
insulator in which < t/U ^ 1, red curve. For < t/U ^ 1, L must lie strictly above 
the U = oo limit and hence L > 2x away from the atomic limit, (b) Hopping processes 
mediated by the t/U terms in the expansion of the projected transformed operators 
in terms of the bare electron operators (see Eq. (flB]) ). As a result of the t/U terms 
in Eq. (|16p . the low-energy theory in terms of the bare fermions does not preserve 
double occupancy. The process shown here illustrates that mixing between the high 
and low-energy scales obtains only if double occupancy neighbours a hole. In the exact 
low-energy theory, such processes are mediated by the new degree of freedom, ipi, the 
charge 2e bosonic field which binds a hole and produces a new charge e excitation, the 
collective excitation in a doped Mott insulator. 

in the atomic limit. In strong coupling, finite hopping with matrix element t creates 
double occupancy, and as a result empty sites with weight t/U. Such empty sites with 
fractional weight contribute to L as shown first by Harris and Lange[l]. In fact, every 
order in perturbation theory contributes to L. Consequently, when < t/U ^ 1, L is 
strictly larger than 2x. Such hopping processes or quantum fluctuations do not affect 
the number of electrons that can be added to the system, however. Since Uh remains 
fixed at 2x, in a real doped Mott system, L/uh > 1. Consequently, in contrast to a Fermi 
liquid, simply counting the number of electrons that can be added does not exhaust the 
available phase space to add a particle at low energy. That is, addition states that do 
not have the quantum numbers of an electron must exist as illustrated in Fig. ([6]). Since 
the number of ways of adding a particle exceeds the number of electrons that can be 
added, the additional states must be gapped to the addition of an electron. This gap can 
manifest itself straightforwardly as a depression in the density of states at the chemical 
potential or more subtely as a reconstruction [5H] of the non-interacting Fermi surface, 
for example, one that has electron (or hole) pockets that shrink in size as the doping 
decreases. Numerical simulations show[59j that such reconstructions do not necessitate 
broken symmetry but obtain entirely from the strong correlations in a doped Mott 
insulator. In either case, the one-to-one correspondence between the excitation spectrum 
in the free and interacting systems breaks down. In doped Mott systems, this breakdown 
arises entirely from spectral weight transfer. While it has been known for some time[l] 
that L > 2x aX strong coupling in a doped Mott insulator based on the Hubbard model, 
that this simple fact implies a pseudogap (whose dimensional dependence is discussed 
in Section 3.2) has not been deduced previously. Thus, additional degrees of freedom at 
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Figure 6. Schematic depiction of electronic models based on the magnitude of L/rih- 
In a Fermi liquid, L/uh — 1. L/rih > 1 necessarily leads to a break down of the Fermi 
liquid picture as new degrees of freedom not made out of the electrons are needed. 
L/uh > 1 appears to be the generic way in which Fermi liquid theory breaks down in 
a doped Mott insulator. 



low-energy, not made out of the elemental excitations, emerge in a low-energy reduction 
of a doped Mott insulator at strong coupling. Note, it is only the spectral weight in 
excess of 2x that creates the new physics. Hence, the physics governed by the new 
degrees of freedom has nothing to do with gauge [601 IGTI [62] fields that engineer the 
no double occupancy constraint in standard treatments of doped Mott insulators. The 
physics referred to here is precisely the part thrown away in such treatments. We 
refer to this contribution as the dynamical part of the spectral weight. This dynamical 
part of L arises through exchanges with doubly occupied sites. Hence, although t/U is 
approximately 1/10 in the cuprates, the t/U corrections to L must be retained as they 
mediate fundamentally new physics that is quantifiable. Namely, it is through these 
corrections that a Fermi liquid description breaks down as depicted in Fig. (El). As a 
result, one might imagine that the new physics is associated with some new collective 
excitation with charge 2e. The new theory we construct, in which we integrate out 
exactly the high energy scale, has such a degree of freedom which does in fact mediate 
the dynamical part of L. In light of Polchinski's argument |42j . any non- Fermi liquid 
behaviour must emerge from the new collective charge 2e excitation. We show that this 
is in fact the case. 

2.3. Perturbative Approaches 

It is possible to account for the dynamical contribution to the low-energy spectral weight 
using degenerate perturbation theory. While this method does not shed any light on 
the missing degree of freedom at low energy, it does serve to illustrate that the limits 
of f/ — >^ CX3 and the thermodynamic limit, oo, do not commute. In fact, it is this 

lack of commutativity that gives rise to L > 2x. We review this method here as it does 
serve to motivate our eventual analysis. 

The goal of perturbative approaches [63 l [64 l [65l [66l [67 l [68] in this context is to bring 
the Hubbard model into diagonal form with respect to double occupancy. As with any 
matrix diagonalization problem, the new basis which makes double occupancy a good 



Mottness in High- Temperature Copper- Oxide Superconductors 



12 



quantum number involves some linear combination of the old states. The subtlety that 
this introduces is that the no double-occupancy condition applies only to the transformed 
fermions not to the original bare electrons. This is an oft-overlooked fact that has led to 
much confusion over what precisely the accepted low-energy reduction of the Hubbard 
model, namely the t-J model[69], entails. We review the derivation with an eye on 
isolating the processes which lead to dynamical spectral weight transfer. Let fi„ be 
the dressed operators which make the Hubbard model block diagonal with respect to 
"double occupancy." Following Eskes et al. [63j, for any operator O, we define O such 
that O = 0(c) and O = 0(/), simply by replacing the Fermi operators Cjo- with the 
transformed fermions fia- To block diagonalise the Hubbard model, 

H[f]^e'^f^H[f]e-'^f\ (10) 

one constructs a similarity transformation S[f] which connects sectors that differ by at 
most one 'fictive' doubly occupied site, that is, a doubly occupied site in the transformed 
basis. To lowest order, 



5« = i - T_i) . (11) 



where 



T+i = -^Yl dij^iaflfjail - rij^), (12) 

which increases the quantum number V = J2i^it'^ii by one. Likewise, T_i = (T_|_i)^ 
decreases V by one. In the new basis, [H, V] = 0, implying that double occupation of 
the transformed fermions is a good quantum number, and all of the eigenstates can be 
indexed as such. This does not mean that [H, V] = 0. If it were, there would have 
been no reason to do the similarity transformation in the first place. V, and not V, 
is conserved. Assuming that V is the conserved quantity results in a spurious local 
SU(2)[70l [7l] symmetry in the strong-coupling limit at half-filling. 

To expose the dynamical contribution to spectral weight transfer, we focus on 
the relationship between the physical and transformed fermions. As expected in any 
degenerate perturbation scheme, the bare fermions, 

Cj(T 6 /jo-C — fia Yf ^ ^ [(^i<? ^ia^fja fjafiafia ~l~ fiaficrfja i (13) 

are linear combinations of the multiparticle states in the transformed basis. We invert 
this relationship to find that 

fia — Cicr + ~ y ] Qij-^ija (14) 



where 



Y 



lUi^ - n. 



Since the low energy theory is captured by the sector in the transformed basis which 
has no double occupancy, it is most relevant to focus on the form of the projected 
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transformed fermions. Using the relations above, we find that as expected, the projected 
transformed fermions 

j 



+ (1 - Uj^) (c]^Ci„ - Cj^ct) d^] (16) 



involve double occupancy in the bare fermion basis. Here = —V^ = 1 and 
bi = J2ja VaCiaCja where j is summed over the nearest neighbors of i. The projected 
bare fermion, (1 — nia)cicr, yields the 2x sum rule, whereas it is the admixture with 
the doubly occupied sector that mediates the t/U corrections. A process mediated 
by these terms is shown in Fig. ([5]). This can be seen more clearly by computing L 
directly using Eq. (fT3l) . The standard treatment [60 l IGT I [62 l [69] of the t-J model ignores 
the dynamical corrections as a hard projection scheme is implemented in which the no 
double occupancy condition applies not only to the transformed but also to the bare 
fermions. As we have pointed out in the introduction, the physics left out by projecting 
out double occupancy is important because it tells us immediately that L/uh > 1 as 
can be seen from the expression for L: 

L^2((l-n,0(l-n,J) (17) 

2t 

= 2((1 - ni|)(l - nil)) + 77 51 9ij{fl [{nja - nia)fj„fj^fi„fia 

+ f\-ah.ir)\ (1 - n,^) + /i-c.) (18) 

As is evident, 2((1 — ?T.j|)(l — ?^jj,)) = 2x in the projected Hilbert space of the dressed 
fermions which corresponds to the 2x sum rule of the static part in the low energy 
spectral weight. The dynamical part of L arises from the tjV corrections. It is 
these corrections that prevent the operator in Eq. flTB]) from being regarded as a free 
excitation. Rather it describes a non- Fermi liquid {L/nh > 1). 

However, it is particularly cumbersome to extract physical insight from the 
canonical transformation method. The primary reason is that any information regarding 
the bare fermions requires that the similarity transformation be undone when any 
experimentally relevant quantity is calculated. Consider for example the electron 
spectral function. In the hard projected version of the t — J model[60l [611 [62l [69] . 
the electron spectral function is assumed to be given by simply the time-ordered 
anticommutator of the transformed fermions. However, Eq. f[T3|) illustrates that this is 
not so. In actuality, the single-particle Green function, 

G(k,^) = -iFT{TciJ^,) = -iFT{ThJl) 

-f iFT—^gik{T {hka — ni^) fka f la f tafia + fiafiafka fL) 

+ iFT—^gjkiTfia inka-njs.)fkaflafjafja + fjafjafka ), 



(19) 
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has t/U corrections in the transformed basis. Although these corrections are naively 
down by a factor of t/U relative to the projected part, their contribution cannot 
be ignored because it is in these corrections that the explicit non-Fermi liquid 
behaviour is hidden. To illustrate, all calculations on the t — J model of the single 
hole problem [721 1731 17^ 175] in a quantum antiferromagnet yield a finite value of Z 
proportional to J/t. However, in the Hubbard model, the situation is not as clear. 
Simulations on finite-size systems reveal that Z oc L^" where a > and hence tends to 
zero as the system size increases. While this calculaton is not conclusive, it is consistent 
with the fact that similar dynamical mean-field treatments of the t-J and Hubbard 
models at finite doping yield a finite conductivity as T ^ in the t — J model[76l|77j but 
a vanishing value in the Hubbard model as T 0[78]. These differences are summarized 
in Table ([1]). The most striking results in Table ([T]) are those for the exponents governing 
the asymptotic fall-off of the density correlations, ac, and momentum distribution 
functions (9) in the t — J (with t = J, the supersymmetric point) and Hubbard models 
in d = 1. Here these quantities can be obtained exactly [791 EDI [HI] for both models 
using Bethe ansatz. In the d=l Hubbard model, the exponent 6 approaches |80[ [8T] 
1/8 asympotically as f/ — oo for any filling. By contrast in the t — J model[79], 6 is 
strongly dependent on doping with a value of 1/8 at half- filling and vanishing to zero as n 
decreases. More surprising, the exponent etc remains pinned [79l [80] at 2 for the [/ — >■ oo 
limit of the Hubbard model at any filling. In fact, at any value of U, ac = 2[79l [ED] 
in the dilute regime of the Hubbard model in = 1. In the t — J model [79] {t = J), 
ac starts at 2 at n = 1 and approaches a value of 4 at n = 0. Note, Oc = 4 is the 
non-interacting value. That is, in d = 1 in the dilute regime, density correlations decay 
as in the U —>■ oo Hubbard model and as r""^ in the t — J (t=J). This discrepancy 
is a clear indicator that relevant low-energy physics is lost if double occupancy of bare 
electrons is projected out in the parameter range considered here. Supposedly, this is 
captured by the t/U corrections to the electron operators in Eq. ( fT6l) . In the procedure 
we outline in the Section HI, we show that all of the physics described by the string of 
operators in the t/U corrections (Eq. (fT6l) is described by a single charge 2e bosonic 
field. 

3. New Theory: Hidden Charge 2e Boson 

The Wilsonian program for constructing a low-energy theory is to explicitly integrate 
over the fields at high energy. The theory that results from this procedure should contain 
all the physics at low energy. In the context of the Hubbard model, it should explicitly 
tell us that L/uh > 1, a key defining feature of a gapped phase without symmetry 
breaking. We now show how this comes about. 

For the Hubbard model, one has to cleanly associate the physics on the energy 
scale U with an elemental field that can be integrated out either by using fermionic 
or bosonic path-integral methods. The Hubbard model in its traditional form does 
not admit such a treatment. However, one can bring the Hubbard model into the 
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Table 1. Comparison between Hubbard and t-J models. The exponents ac and 9 were 
computed using Bethe ansatz for the supersymmetric {t — J) t—J and Hubbard models 
in d = 1. Respectively, these exponents govern the asymptotic form of the density 
correlations and the momentum distribution functions. Zih denotes the quasiparticle 
weight for a single hole in Mott insulator described by either the Hubbard or t-J 
models, a{T = 0,n > 0.9) the conductivity for fillings exceeding n = 0.9, and Mid-IR 
denotes the mid-infrared band in the optical conductivity 



t — J Hubbard 



ac{n = 0) 


4[7S] 


2180] 





(tte - 4)Vl6a,[79] 


1/8 f7^(X)[80] 


Zih 


finitelZlElEl] 


L-«[75] 


Mid-IR 


none [Ha [83] 


yes [Si [85] 


Luttinger surface at n=l 


none[86l[8l[[89] 


yes[87|[59l[89] 



appropriate form by introducing an elemental field that describes the excitations far 
away from the chemical potential. For hole-doping, this constitutes the excitations in 
the upper-Hubbard band. In our construction, we will extend the Hilbert space of the 
Hubbard model to include an extra degree of freedom which will represent the upper 
Hubbard band. By a constraint, the new field will represent the creation of double 
occupancy. This field will enter the Lagrangian with a mass of U. When the constraint 
is solved, we recover the Hubbard model. However, since the new field has a mass of 
U, the exact low-energy theory is obtained by integrating over this field rather than by 
solving the constraint. Consequently, the new low-energy theory will contain an extra 
degree of freedom having to do with the coupling with the high energy scale. Since the 
constraint field has to do with double occupancy, it must have charge 2e. 



3.1. Bohm/Pines Redux 

There is a great similarity between our treatment of the new collective mode in the 
Hubbard model and the Bohm-Pines derivation of plasmons. Here we briefly introduce 
the approach used by Bohm and Pines[90] to describe the collective excitation of the 
interacting electron gas. Shankar and Murthy[9T] also used a similar approach in their 
dipole analysis of the u = 1/2 quantum Hall state. The basic idea is to re-express 
the Hamiltonian in such a way that the long-range part of the Coulomb interactions 
between the electron is described in terms of collective fields (plasma mode) by enlarging 
the original Hilbert space. After we remove the unphysical states by a constraint, the 
resultant Hamiltonian will transform to an interacting electron gas with only short-range 
Coulomb interactions coupled to a plasma oscillating mode. 

The starting point is the general interacting electron Hamiltonian, 



H= y ^ + 27re2V^— 2 W V -L 

^ 2m ^ 

kij k 

where n is the total electron density, the first term corresponds to the kinetic energy of 
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the electrons, the second term is their Coulomb interaction and the third term is the 
self-energy which represents a uniform positive charge background. 

The collective mode can be described by first enlarging the Hilbert space of the 
original electron gas to include a new set of canonical variables, (tt/c, 6k) such that 
[^fcjTTfc'] = ih6k,k'- The original Hamiltonian becomes, 

ik k 

^^Y.^k-Me'^^^''-^ (20) 

ikl 

Here, is the unit vector along the k direction. The relevant equation can be derived by 
rewriting the Hamiltonian as a non-interacting electron system coupled to an external 
electric field, E(x), 

«=I:5^^(p<-^■4.(x))^^E(x)^ (21) 

such that 

E(x) = v^5]7r_,efce^'^-. (22) 

k 

The corresponding longitudinal vector potential A(x) is 

A(x) = v^EW'" (23) 

k 

Here, both A(x) and E(x) are real and the unphysical states can be removed by the 
constraint, 

nk = nk-i[-^^ ^e-^'^-^' = OVk, (24) 

which was obtained by equating the electric field energy, E(x)^/87r with the electron- 
electron interaction energy. For the last term in Eq. fl20|) . the dominant part comes 
from k = — 1. By defining the plasma frequency 

ul = , 25 

we are able to simplify the Hamiltonian 

2"-^ i: (p) + ^ E • (P. - m<^- 

I k ik 

- ^ E {^k-n^k + ^pke-k) , (26) 
^ k 

which describes the non-interacting electron gas coupled with the plasma mode of 
frequency Up. Here, we have simplified the derivation by assuming the collective modes 
can oscillate with any frequency. In a realistic system, a maximum cutoff frequency. 



Mottness in High- Temperature Copper- Oxide Superconductors 



17 



|0,0;0) 

1,0;0) |0,0;1) |0,1;0) 




|l,0;l) |l,l;0) |0,1;1) 



(0,0) (0,<7) (0,D) 
(a,0) {a -a) (a,D) 
(D,0) {D,cj) (D,D) 



Figure 7. a) Extended Hilbert space which allows a clean coarse-graining on the 
energy scale U. b) Hopping processes which are included in the Lagrangian in the 
extended Hilbert space. 



kc, determined by the electron density, arises so that only the long-range electron- 
electron interaction can be transformed into the plasma mode, and the electron gas 
retains a short-range Coulomb interaction. The magnitude of can be determined self- 
consistently by minimizing the total energy. To summarize, we have mapped the original 
electron-electron interacting Hamiltonian to a non-interacting electron gas coupled to 
the plasma mode. The key trick that made this happen was enlarging the original 
Hilbert space. 



3.2. Charge 2e Boson 

The essence of the Bohm-Pines|90] derivation is that plasmons, as an independent degree 
of freedom, are only apparent when the constraint is relaxed in the extended Hilbert 
space. As we will see, the same is true for a doped Mott insulator. To this end[9l|in], we 
extend the Hilbert space ®j (JF^ ® JF^ J-'d), where J-' denotes a fermionic Fock space. 
In the left frame of Fig. (JTj), we indicate the states of the extended Hilbert space for 
a single site. The new extended states are shown in red. Such states will be removed 
by a suitably chosen constraint which will associate D] with the creation of double- 
occupation. In order to limit the Hilbert space to single occupation in the D sector, we 
will take D to be fermionic. The field D will enter the theory as an elemental field with 
a large (order U) quadratic term and precise interactions with the electronic degrees of 
freedom; the low-energy (IR) theory is obtained by integrating out D. Because D is 
fermionic, a trick is required to associate it with the creation of double occupancy which 
clearly transforms as a boson. Essentially, we will have to fermionize double occupancy. 
This can be done by imposing a constraint on the Di field such that 

A ~ e^^Cii (27) 
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where ^ is a Grassmann field which is normahzed as 

d'^e ee = i. (28) 

The constraint will be imposed through the use of a 5— function and hence will enter 
the action upon a subsequent exponentiation as in the implementation of the constraint 
in the non-linear cr— model. While there is some similarity between 6 and a super- 
coordinate, this association is strictly formal in our case as any dependence on the 
Grassmann parameter 6 disappears. 

For hole-doping, it is the upper Hubbard band that must be integrated out. The 
appropriate hopping processes that must be included in the Lagrangian are given in the 
right-hand panel of Fig. ([7]). The Euclidean Lagrangian in the extended Hilbert space 
which describes the hopping processes detailed above can be written 



c= d^e 



i,a i 

+ UY^ D]Dj - t Qij Cij^^cl^Cj^^ + Djc]^^Ci^^Dj 



+ (L'J^Q,<,Kc„-. + h.c.)\ + H,on\ . (29) 

Here, gij selects out nearest neighbours (note that if we wanted to include next-to- 
nearest neighbour interactions, we need just modify the matrix accordingly), the 
parameter has values V| = 1, = —1, and simply ensures that D couples to the 
spin singlet and the operator Cij^o- is of the form Cj^ o- = 09aij^„ = 99{l — ni _f^){l — nj^^f^) 
with number operators Ui^a = cl^Ci,^. Note that the dynamical terms that appear in the 
Lagrangian exclude, as they must, those sites already singly occupied, as they describe 
the dynamics in the LHB. The constraint Hamiltonian ifcon is taken to be 

i^con = seY,^]{D, - ^9, tc,, J + /i.e., (30) 

3 

where is a complex charge 2e bosonic field which enters the theory as a Lagrange 
multiplier. The constant s has been inserted to carry the units of energy. It could be 
absorbed into the definition of the constraint field There is a natural connection 
between yjj in our theory and a in the non-linear cr-model. Both start their lives as 
Lagrange multipliers but both end up affecting the low-energy physics. At this point, 
there is some ambiguity in the normalization of <y9, but we expect that this will be set 
dynamically. We will find that if a true infrared limit exists, then s must be of order 
the hopping matrix element t. 

Now, as remarked previously, we have chosen the Lagrangian ( l29l) so that this 
theory is equivalent to the Hubbard model. To demonstrate this, we first show that 
once the constraint is solved, we obtain the Hubbard model. Hence, the Lagrangian we 
have formulated is the Hubbard model written in a non-traditional form - in a sense, 
we have inserted unity into the Hubbard model path integral in a rather complicated 
fashion. To this end, we compute the partition function 

Z = J [Do Vc^ VD VD^ Vip Vip^] exp^/o"^'^* . (31) 



Mottness in High- Temperature Copper- Oxide Superconductors 



19 



with C given by (12^ . We note that ip is a Lagrange multipher. As shown previously |10j. 
in the Euchdean signature, the fluctuations of the real and imaginary parts of must 
be integrated along the imaginary axis for (pi to be regarded as a Lagrangian multiplier. 
The (fi integrations (over the real and imaginary parts) are precisely a representation of 
(a series of) 5-functions of the form, 



6 / deD, 



dO Oci tc. 



(32) 



If we wish to recover the Hubbard model, we need only to integrate over Dj, which 
is straightforward because of the 5-functions. The dynamical terms yield 



d^e ee ^(i - ni,_^)4^Q,^ + ^ 4i4^Ci,tCi,i 

i 

j d^e 9eJ2 [(1 - rii_^)cl^Ci^^ + ni_^cl„Ci^^ 

'd^eeej^clc,,^. 

Likewise the term proportional to yields 

= j d^O 66 gijUj^^^cl^^Ci^a + h.c. 
Finally, the two D field terms give rise to 

d'^6 e6j29ij [cl,i4T(4TCi.T + 4i^^i)^i.TCj-,i 

id 

- J d^e 6eJ2gijnj^_^ni_^cl^Cj^a- 



(33) 



h.c. 



(34) 



(35) 



Eqs. flMl) and fl35l) add to the constrained hopping term in the Lagrangian (the 
term proportional to Cjj,o-) to yield the standard kinetic energy term in the Hubbard 
model. Finally, the D^D term generates the on-site repulsion of the Hubbard model. 
Consequently, by integrating over yjj followed by an integration over Di, we recover the 
Lagrangian, 

d^9 99C}iuhh = X! 4,0-^4,0- + -f^^Hubb, 



(36) 



of the Hubbard model. This constitutes the ultra-violet (UV) limit of our theory. In 
this limit, it is clear that the Grassmann variables amount to an insertion of unity 
and hence play no role. Further, in this limit the extended Hilbert space contracts, 
unphysical states such as |1, 0, 1), |0, 1, 1), |1, 1, 1) are set to zero, and we identify |1, 1, 0) 
with |0,0, 1). Note there is no contradiction between treating D as fermionic and the 
constraint in Eq. fl30p . The constraint never governs the commutation relation for D. 
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The value of D is determined by Eq. (1301) only when (f is integrated over. This is 
followed immediately by an integration over D at which point D is eliminated from the 
theory. 

The advantage of our starting Lagrangian over the traditional writing of the 
Hubbard model is that we are able to coarse grain the system cleanly for U ^ t. 
The energy scale associated with D is the large on-site energy U. Hence, it makes sense, 
instead of solving the constraint, to integrate out D. The resultant theory will contain 
explicitly the bosonic field, (f. As a result of this field, double occupancy will remain, 
though the energy cost will be shifted from the UV to the infrared (IR). Because the 
theory is Gaussian, the integration over Di can be done exactly. This is the ultimate 
utility of the expansion of the Hilbert space - we have isolated the high energy physics 
into this Gaussian field. As a result of the dynamical term in the action, integration 
over D will yield a theory that is frequency dependent. The frequency will enter in the 
combination u U which will appear in denominators. Since U is the largest energy 
scale, we expand in powers oico/U; the leading term yields the proper a; = low-energy 
theory. Since the theory is Gaussian, it suffices to complete the square in the D-Held. 
To accomplish this, we define the matrix 

and hi = Y^jbij = Y.j,a 9ij(^j,o-VaCi-cT- At zero frequency the Hamiltonian is 



where 



t^^,,,. .... 



- s E vh^c.^i + ^ E ^\{M-'),,b, + h.c. (38) 

j i,j 

which constitutes the true (IR) limit as long as the energy scale s is not of order U. As 
we have retained all powers of t/U, Eq. ( 138|) is exact. If s ~ 0{U) then we should also 
integrate out ^ this integration is again Gaussian and can be done exactly; one can 
easily check that this leads precisely back to the UV theory, the Hubbard model. Hence, 
the only way in which a low-energy theory of the Hubbard model exists is if the energy 
scale for the dynamics that cp mediates is 0{t). This observation is significant because it 
lays plain the principal condition for the existence of an IR limit of the Hubbard model. 
Since the order of integrations we have performed here does not matter, integration 
over in the path integral for the action corresponding to Eq. (!38ll also yields the 
Hubbard model. As we have shown elsewhere [10] the Wick rotation must be taken 
into consideration to complete the Gaussian integration over cpi. Finally, as we have 
shown previously[TO], the theory derived here could easily have been constructed in 
terms of the Hubbard operators, ^ and r]. This offers no further complication. In so 
doing, the spin-spin interaction which arises from the first term in Eq. ( l38l) would have 
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been projected onto a subspace which prohibits double occupancy. That is, it would 
be equivalent to the spin-spin interaction in the standard t — J model. Since double 
occupancy still survives at low energies and is mediated by the ipi terms, such a rewriting 
of the low-energy Hamiltonian is strictly optional. 

To fix the constant s, we determine how the electron operator transforms in the 
exact theory. As is standard, we add a source term to the starting Lagrangian which 
generates the canonical electron operator when the constraint is solved. For hole-doping, 
the appropriate transformation that yields the canonical electron operator in the UV is 



+ h.c. 



However, in the IR in which we only integrate over the heavy degree of freedom, Di, 
the electron creation operator becomes 

t s 

^ (1 - n^,-a)cl^ + V^-^hCi^^^ + K^'/^Iq,-^ (39) 

to linear order mt/U. This equation resembles the transformed electron operator in Eq. 
f ITB]) . In fact, the first two terms are identical. The last term in Eq. flTB]) is associated 
with double occupation. In Eq. (139|) . this role is played by ifi. Demanding that Eqs. 
( [T6|l and (139!) agree requires that s = t, thereby eliminating any ambiguity associated 
with the constraint field. Consequently, the complicated interactions appearing in Eq. 
( IT6|l as a result of the inequivalence between Cia and fia are replaced by a single charge 2e 
bosonic field (fi which generates dynamical spectral weight transfer across the Mott gap. 
The interaction in Fig. ([5]), corresponding to the second-order process in the term flbi, 
is the key physical process that enters the dynamics at low-energy. That the dynamical 
spectral weight transfer can be captured by a charge 2e bosonic degree of freedom is 
the key outcome of the exact integration of the high-energy scale. This bosonic field 
represents a collective excitation of the upper and lower Hubbard bands. Hence, we have 
successfully unearthed the extra degree of freedom associated with L/uh > 1 physics in 
a doped Mott insulator. 



3.3. Half-Filling: Mott gap and antiferromagnetism 

There is an important simplification that obtains at half-filling that points to one of 
the potential uses of this theory: the dynamical mode that generates the Mott gap. 
Recall our ultimate task was to integrate out the degrees of freedom far away from the 
chemical potential. At half-filling, both the lower and upper Hubbard bands are at 
high energy. Hence, both must be integrated out to obtain the true low-energy theory. 
However, at present, we have only integrated out the UHB. Because the integration of 
the UHB is not equivalent to simply integrating out double occupancy, a trivial particle- 
hole transformation does not help us to determine how the low-energy theory should be 
formulated in this case. 

Ultimately we have to introduce two new fermionic fields D and D associated 
with the double occupancy and double holes, respectively. To proceed, we consider the 
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Lagrangian, 

/:hf^= Jd^e iD^i) + iDD -^{D^D - DD^) 
+ ^D^b + ^ebD + h.c. + sOif^D - 0C|q) 

+ SOif^D - ec\cl) + h.c] , (40) 

contains the two constraint charge ±2e bosonic fields, (pi (charge 2e) and (pj (charge 
— 2e) which enter the theory as Lagrange multipliers for the creation of double occupancy 
and double holes, respectively. Similar to the previous result, if we first integrate out 
both the bosonic fields ipi, (pi and then Z)j, Dj, the Hubbard model is obtained and 
the generalised theory EqJlO] yields the correct UV limit. In deriving the UV limit, it 
is crucial that the operators representing the creation of double occupancy and double 
holes remain in the order indicated in the Lagrangian. However, a different IR limit is 
obtained if we first integrate out Di and Dj, 



2 

sip^ — s*(pj c-fCj + h.c, (41) 

where L± = i-^ ± ^. This theory is exact and hence should contain the source of the 
Mott gap. This Lagrangian is invariant under the transformation Ci^a exp(zQ-Rj)c|o., 
y^i ^ 'Pi and s ^ s. This invariance refiects the symmetry between the double 
occupancy and the double hole in the system at half filling. In contrast to the doped 
case as in Eq. (l38l) . no M. matrices appear in the IR theory at half filling. Consequently, 
we arrive at a closed form for the low-energy theory at half-filling in which no bare field 
has dynamics. The 6^6 + terms include a spin-spin interaction as well as a three-site 
hopping term. However, at half-filling, the three-site hopping term vanishes. 

It is interesting to note that Eq. fHT]) . as an exact low-energy theory of the 
Hubbard model, is not equivalent to the Heisenberg model. Only the b'^h terms yield 
the Heisenberg model. Hence, we consider the separation, 

-^Mott = + U^^LZ^h + -fh'^L-^h (42) 

which explicitly removes the spin-spin term from the low-energy Lagrangian. As will 
be seen, the dynamics leading to the Mott gap can be constructed entirely from Z^Mott- 
That C}^ is not equivalent to the Heisenberg model is not surprising for three related 
reasons. First, explicit evaluation of the path integral over the new degrees of freedom 
must regenerate the original Hubbard model. Hence, there must be something left over 
once we subtract the Heisenberg terms from the low-energy theory of the Hubbard 
model. Second, as pointed out previously [70l ITT] , the Heisenberg model has a local 
SU(2) symmetry not the global SU(2) symmetry of the Hubbard model. Hence, the 
true low-energy theory of the Hubbard model at half-filling must have other terms that 
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break the local SU(2) and reinstate the global SU(2) symmetry. As we have shown 
previously jlO], all of the terms involving the cpi and ipi degrees of freedom restore the 
global SU(2) symmetry. Consequently, the emergence of the new local SU(2) symmetry 
is a function entirely of projection onto the singly-occupied electron subspace. Third, at 
half-filling, the Hubbard model possesses a surface of zeros [86l EZ] of the single-particle 
Green function along a connected surface in momentum space, whereas the Heisenberg 
model does not. The absence of such a surface of zeros, the Luttinger surface, is also 
a function of projection. In fact, all of these three failures of the Heisenberg model 
arise from hard projection, which Eq. (fT6|) shows is not correct even at half-filling. The 
non-trivial implication of the zero surface is that the real part of the Green function, 

RAK 0) = - de'M^ - r de^^^ (43) 



e JA+ e 

vanishes. Here Afj{k,e) is the single-particle spectral function which we are assuming 
to have a gap of width 2A symmetrically located about the chemical potential at e = 0. 
Because A(k, e) > away from the gap, and e changes sign above and below the gap, 
Eq. fH3|) can pass through zero. For this state of affairs to obtain, the pieces of the 
integral below and above the gap must be retained. Projected models which throw away 
the UHB fail to recover the zero surface. 



3.3.1. Mott Gap What Eq. ( 1411) makes clear is that all the information regarding 
the surface of zeros is now encoded into the bosonic fields yjj and (^j. While numerical 
met hods [921193] exist which lead to a Mott gap, there has been no explicit demonstration 
of the dynamical degrees of freedom that ultimately generate this gap. The bosonic 
degrees of freedom in Eq. (HT!) solve this problem. 

To proceed, we transform to frequency and momentum space and focus on a 
square lattice as in the cuprates. Defining ip{t) = J doj e~'^'^^{p^, the energy dispersion, 

_ 4 cos(A;^a/2) cos(p^a), and the Fourier transform of 6j, 

= Y.^?^ Ck/2+p,TCk/2-p,i, (44) 

p 

we arrive at the exact working expression, 

/•hf I I I I I 

= " (c.-^/2)^-''^^"''^ + (Z7TW2) ^-"''^^""•'^ 

+ K'H^)^i.k + 5*«?H^)^-.,k) 

X (ck/2+p,TCk/2-p,i)c. + ^-c- (45) 
for the low-energy Lagrangian where we have suppressed the implied integration over 
frequency and introduced the coupling constants. 
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which play a central role in this theory. They, in fact, will determine the spectral weight 
in the lower and upper Hubbard bands, respectively. 

Although Eq. ( H5i) lacks any kinetic terms, an analysis similar to that used by 
Polchinski [42j in the context of Fermi liquids can be used. The key point in that 
argument is that all renormalizations are towards the surface in momentum space where 
the spectral weight resides. In a theory of weakly coupled constituents, setting the 
coefficient of the quadratic terms in a Lagrangian would determine their dispersion. As 
is evident from Eq. ( l45l) . the coupling constants for all of the quadratic terms never 
vanish. All else being equal, this implies that there is no spectral weight anywhere. What 
we have shown [95] is that the terms in which the bosonic degrees of freedom and the 
fermions are coupled determine where spectral weight resides. In fact, it is these terms 
that should properly be regarded as quadratic. To this end, we make the transformation, 
iPuj \Jl - 2u/U (p^ + 2ujIU (p^, and (cc)^; \Jl- 4co'2/f/2 (^^^.^^ which 

recasts the Lagrangian as 

li i9 i9 ^ 



+ S7j'H^)</'l,kCk/2+p,c^/2W,TCk/2-p,^/2-..'a (48) 
+ 5*7^''^ ('^)<^-tj,kCk/2+p,w/2+a;',TCk/2-p,cj/2-a;',i 

+ h.c. (49) 



where 



. 1±^^l±^,JVZ^j. (50) 

The vanishing of the coefficients 7p and 7p determine where the spectral weight lies. 
Consider initially k = so that the dispersion simplifies to 4 cosp^. When oj = U /2 
{oj = —U/2), J2fi cospf^ = defines the momentum surface along which 7p (7p) vanishes. 
This surface corresponds to the diamond ap = {ap^, ±n — apx) as depicted in Fig. ([9]). 
These features define the center of the LHB {—U/2) and UHB {U /2). At each momentum 
in the first Brillouin zone (FEZ), spectral weight develops at two distinct energies. This 
state of affairs obtains because 7p = between U /2 — At < u <U /2 + At and 7p = for 
—U /2~ At <uj < —U /2 + At for each momentum state in the FEZ. Within each energy 
range, the associated operator, which is of the form Lp'^cc (UHE) or (pec (LHE), can be 
viewed as a quadratic kinetic term. The (0, 0) point corresponds to the lowest energy 
state in each LHE and UHE, that is, uj = —U /2 — 4t and u = U/2 — At whereas (tt, tt) 
sits at the top of each band at = —U/2 + At and u = U/2 -\- At. Consequently, at 
each momentum, the splitting between the turn-on of the spectral weight in the UHE 
and LHE is U. When U > W, each momentum state lacks spectral weight over a 
common range of energies. As a consequence, a hard gap opens in the spectrum. This 
is the Mott gap (for the composite excitations not the electrons), and its origin is the 
emergence of composite excitations described by the operators ip'^cc (UHE) and ipcc 
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Figure 8. a) Diamond-shaped surface in momentum space where the particle 
dispersion changes sign. b) Turn-on of the spectral weight in the upper and 
lower Hubbard bands as a function of energy and momentum. In the UHB, the 
spectral density is determined to 7p while for the LHB it is governed by 7p. The 
corresponding operators which describe the turn-on of the spectral weight are the 
composite excitations (^s'l'cc (UHB) and (fcc (LHB). 



(LHB), which we loosely interpret as bound states. As our analysis thus far is exact, 
we conclude that in the absence of any symmetry breaking, the coefficients 7p and 7p 
determine the dispersion for the excitations that comprise the here-to-fore undefined |4 7] 
UHB and LHB. Inclusion of the center-of-mass momentum k simply shifts the value of 
the momentum at which the dispersion changes sign, thereby keeping the Mott gap in 
tact. 

Thus far, we have established the Mott gap in terms of a set of composite excitations 
which are orthogonal in that they never lead to a turn-on of the spectral weight in the 
same range of energies. Ultimately, we would like to know the spectral function in 
terms of the original electron degrees of freedom. The lack of any derivative terms in 
the action with respect to ipi implies that we can treat ip as a spatially homogeneous 
field. A priori, gradient terms with respect to (fi are possible. However, such terms 
are absent from the exact low-energy theory as such terms would indicate the presence 
of a freely propagating bosonic degree of freedom at half-filling. It is precisely because 
such terms are absent that we were able to identify that the only propagating degrees 
of freedom at half-filling are gapped excitations. 

To proceed, we rewrite the coefficient of the boson-fermi terms as 



In this treatment, any non-trivial dynamics leading to the Mott gap will arise only 



A{k,u,(f),(f)) 



s(0t - 0) 




and 




(51) 




Performing the Wick rotation, (p ^ icf) and 
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0* — > i(f)*, we recast the single-particle electron Green function as 

G{k,u) = J d(j)J (i0G(fc,a;,0,0)exp-/^^M°" (52) 

where 

^' = I^^|2^-A - ^^^^ 

|A(A;, uj, (p, + Id 

At first glance, the Green function seems to have a vanishing imaginary part. However, 
because of the i5 in the gap function, A(A;, a;, 0, 0), the imaginary part of the Green 
function 



^G{k, oj, 0, 0) = lim {U - 2ujf + 5^ (f/ + 2u;)2 + 5^ 

5 



X 



+ (2A(0 + 0) + 52) 52 + 0(5 
(f/-2cu)2(?7 + 2cu)2 



5 



5{A) (54) 



is explicitly non-zero. Here 



A 



- Auj^ - 2ak{U + 2uj) 0+ f/^ - 4cu^ - 2afc([/ - 2u;) 



S = 20(2cu + ftfc) + 20(2cj + Ofe) (55) 

To complete the calculation, we performed the cpi and if integrations numerically. The 
results in Fig. (Q clearly show that a Mott gap exists and the spectral weight is 
momentum dependent. At (0,0), the spectral weight lies predominantly in the LHB 
whereas at (vr, vr) it lies in the UHB. Consequently, the real part of the Green function 
must change sign on some momentum surface between these two limits. The location of 
the zero surface is the Fermi surface of the non-interacting system as it must be for the 
half-filled system with particle-hole symmetry. We find then that the Mott gap arises 
from the dynamics of the two charge |2e| bosonic fields. This is the first time the Mott 
gap has been derived dynamically, in particular by a collective degree of freedom of the 
lower and upper Hubbard bands. Relative to the gap in the spectrum for the composite 
excitations that diagonalise the fermion-boson terms in Eq. (l49l) . the gap in the electron 
spectrum is larger. This is not surprising as the bare electrons do not have unit overlap 
with the composite excitations. In addition, the momentum dependence of the spectral 
function is identical to that obtained by dynamical mean-field calculations |92j thereby 
lending crecedence to such cluster [M]. 

An open question that this analysis provokes is whether or not the turn-on of the 
spectral weight in a Mott insulator is governed by a fixed point. If so, then in analogy 
with the Fermi liquid analysis ^2]. all the interactions except those that govern the 
turn-on of the spectral weight should be irrelevant. That is, the (p^cc and the ipcc 
terms represent a natural theory. Indeed, the analysis presented above demonstrates 
that the spin-spin term has nothing to do with the turn-on of the spectral weight, as 
foreshadowed by Mott[T]. Namely, the gap in the spectrum at half-filling is independent 
of ordering. While a naive scaling analysis suggests that the spin-spin interaction is 
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Figure 9. Electron spectral function corresponding to ^Mott for U = 8t. The gap 

here is generated entirely from the dynamics of the charge 2c bosonic fields that emerge 
from integrating out the upper and lower Hubbard bands at half-filling. 



indeed sub dominant, we have been unable to compute the /3-function to show that a 
true fixed point underhes the physics at half-fiUing. Such a computation stands as a 
true challenge for Mottness. 

Nonetheless, antiferromagnetism with an ordering wave-vector of (tt, tt) can also 
be understood within this formalism. Within this theory, there is a natural candidate 
for the antiferromagnetic order parameter, namely Bij — (gijiplci^-^Cj^i). The vacuum 
expectation value of this quantity is clearly non-zero as it is easily obtained from 
a functional derivative of the partition function with respect to 7p. That this is 
the relevant order parameter instead of the traditional one follows from the fact 
that the spin-spin interaction and all higher-order operators contained in |6p are at 
least proportional to (a the lattice constant) and hence are subdominant to the 
composite terms. Hence, a non-traditional order parameter must govern the turn-on 
of antiferromagnetism. We advocate that Bij characterizes the antiferromagnet that 
describes the strong-coupling limit of the Hubbard model and as a consequence the 
insulating state of the cuprates. An antiferromagnet of this kind has no continuity with 
the antiferromagnet at weak coupling because it is mediated by the collective mode (p 
or (f. Hence, both the Mott gap and subsequent antiferromagnetic order emerge from 
composite excitations that have no counterpart in the original UV Lagrangian but only 
become apparent in a proper low-energy theory in which the high-energy degrees of 
freedom are explicitly integrated out. Away from half- filling, a similar state of affairs 
obtains. 

4. Hole Doping: Experimental Consequences 

The charge 2e boson has much to tell us about the normal state of a doped Mott 
insulator. Here we compute the electron spectral function, the specific heat, the thermal 
conductivity, the optical conductivity as well as the dielectric function. In each of these, 
the charge 2e boson produces a distinct signature that accounts for the anomalies of the 
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doped state of a Mott insulator. 
4-1- Spectral Function: Pseudogap 

Since we have demonstrated that Z^Mott captures the strong-couphng physics of the Mott 
insulating state, we focus on the evolution of this theory with doping. The lack of any 
gradient terms in the action with respect to the bosonic fields and the absence of any 
bare dynamics associated with ipi suggests that we can treat ipi as a homogeneous field. 
Further, since we are not interested in the dynamics on the Mott scale, we treat 9? as a 
static field. Its sole role is to mix the subsectors which differ in the number of doubly 
occupied sites. Consequently, our results are valid provided that uj < U and U ^ t. 
Under these assumptions, the single-particle electron Green function 

G{k,uj) = -iFT{Tc^{t)c]{0)), (56) 

can be calculated rigorously in the path-integral formalism as 

Gik,uj) = -iFT I [DcpfllDcp,] I [Dc*][DQ]Q(t)c*(0)exp-/^[=''^l'^*, (57) 

where FT refers to the Fourier transform and T is the time-ordering operation. The 
explicit spin-spin term is not contained in Z^Mott- This term will also be dropped in the 
doped case because even in this limit, the spin-spin term is subdominant (in a naive 
continuum limit sense) to the other interactions in C This state of affairs arises because 
spin-spin term in 6^6 contains four spatial derivatives, whereas the ip^b term contains 
only two. As a result, all of the physics we present below is associated with the charge 
rather than the spin degrees of freedom. The continuity of the analysis with that of the 
half-filled system raises the question that perhaps a fixed point at half-filling persists 
to finite doping as well in which only the fermion-boson terms are relevant. While our 
analysis is highly suggestive that such a state of affairs obtains, the possible existence 
of such a fixed point remains a conjecture as of this writing. 

To proceed, we will organize the calculation of G{k,uj) by first integrating out the 
fermions (holding fixed) 

Gik.u) = I [Dip*][Dif]FT [Z}c*][Z}Q]Q(t)c*(0)exp-/^[^''^]'^*) (58) 

where now 

ia \ ^ / ka 

2t 

+ - -^)C-feiCfct + C.C. (59) 

k 

The effective Lagrangian can be diagonalized and written in terms of a collection 
of Bogoliubov quasiparticles[TT]. The remaining ip integration can then be done 
numerically to obtain the spectral function. 

The spectral functions in Figs. (fTOj) and (fTTj) exhibit four key features. First, 
there is a low-energy kink in the electron dispersion that is independent of doping. The 
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Figure 10. (a) Spectral function for filling n = 0.9 along the nodal direction. The 
intensity is indicated by the color scheme, (b) Location of the low and high energy 
kinks as indicated by the change in the slope of the electron dispersion, (c) The energy 
bands that give rise to the bifurcation of the electron dispersion. 



low-energy kink occurs at roughly 0.2t ^ lOOmeV. By treating the mass term for the 
boson as a variable parameter, we verified that the low-energy kink is determined by 
the bare mass. In the effective low-energy theory, the bare mass is t'^/U. This mass is 
independent of doping. Experimentally, the low-energy kink|96j does not change with 
doping. Consequently, the charge 2e bosonic field provides a natural mechanism for the 
kink that is distinct from the phonon schemes that have been proposed [96 j. 

Second, a high-energy kink appears at roughly 0.5t ~ 250meV which closely 
resembles the experimental kink at 300mel^[97]. Cluster [98] and exact diagonalization 
methods [99] also find a high-energy kink. At sufficiently high doping (see Figs. ffTTk ) and 
flllb)). the high-energy kink disappears. Third, the electron dispersion bifurcates at the 
second kink. This is precisely the behaviour that is seen experimentally [97] . The energy 
difference between the two branches is maximum at (0, 0) as is seen experimentally. A 
computation of the spectral function at f/ = 20t and n = 0.9 reveals that the dispersion 
as well the bifurcation still persist. Further, the magnitude of the sphtting does not 
change, indicating that the energy scale for the bifurcation and the maximum energy 
splitting are set by t and not U. The origin of the two branches is captured in Fig. (ITOb ). 
The two branches below the chemical potential correspond to the standard band in the 
LHB (filled squares in Fig. ffTOb ) on which cp vanishes and a branch on which 
(filled circles in Fig. (fTOb ). Simulations on the Hubbard model clearly resolve either the 
low-energy feature [92l [93l [9l] or the high-energy kink [981199] . In the studies showing the 
hihg-energy kink, the low-energy feature is not discernible [9 8 1 [99] . What is new here is 
that both features (but with drastically different intensities as is seen experimentally) are 
captured. The two branches indicate that there are two local maxima in the integrand in 
Eq. ( l58l) . a feature not captured by a saddle-point approximation. Above the chemical 
potential only one branch survives. The split electron dispersion below the chemical 
potential is consistent with the composite nature of the electron operator dictated by 
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Figure 11. Spectral function for two different fillings (a) n = 0.8 and (b) n = 0.4 along 
the nodal direction. The absence of a splitting in the electron dispersion at n = 0.4 
indicates the bifurcation ceases beyond a critical doping. The spectral functions for 
two different values of the on-site repulsion, (c)C/ = lOi and (d)C/ = 20i for n — 0.9 
reveals that the high-energy kink and the splitting of the electron dispersion have at 
best a weak dependence on U . This indicates that this physics is set by the energy 
scale t rather than U . 

Eq. ( l39l) . At low energies, the electron is a linear superposition of two states, one the 
standard band in the LHB described by excitations of the form, c|^(l — ni^) and the 
other a composite excitation consisting of a bound hole and the charge 2e boson, Cia^\. 
The former contributes to the static part of the spectral weight transfer (2a;) while the 
new charge e excitation gives rise to the dynamical contribution to the spectral weight 
transfer. Because the new charge e state is strongly dependent on the hopping, it should 
disperse as is evident from Fig. (fTTj) and also confirmed experimentally [97] . 

The formation of the composite excitation, QctV'^ is the new dynamical degree of 
freedom in the doped theory. This dynamical degree of freedom has no counterpart 
in the UV scale. Such a binding of a hole and the charge 2e bosonic field leads to a 
pseudogap at the chemical potential, as evidenced by the absence of spectral weight at 
the chemical potential for both n = 0.9 and n = 0.8. Non-zero spectral weight resides 
at the chemical potential in the heavily overdoped regime, n = 0.4, consistent with the 
vanishing of the pseudogap beyond a critical doping away from half-filling. Because the 
density of states vanishes at the chemical potential, the electrical resistivity diverges as 
T — > 0. Such a divergence is shown in Fig. (fT2h) and is consistent with our previous 
calculations of the dc resistivity using a local dynamical cluster method [78]. In the 
absence of the boson (Fig. (|T2b )). localization ceases. Although this calculation does 
not constitute a proof, it is consistent with localisation induced by the formation of the 
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Figure 12. (a)dc electrical resistivity as a function of temperature for n = 0.9 (b) 
Setting the bosonic degree of freedom to zero kills the divergence of the resistivity as 
T — !■ 0. This suggests that it is the strong binding between between the fermionic and 
bosonic degrees of freedom that ultimately leads to the insulating behaviour in the 
normal state of a doped Mott insulator. 

bound composite excitation, Ciafl- This state of affairs obtains because the boson has 
no bare dynamics. It may acquire dynamics at 0(t^/[/^) as can be seen by expanding 
the Ai matrix in Eq. (l38ll . 

Such bound-state formation lays plain how the strong coupling regime of a doped 
Mott insulator depends on the dimensionality, the doping and the connectivity of the 
lattice. As the charge 2e boson is a local degree of freedom with no bare dynamics, an 
analogy with bound state formation by a local potential is warranted. It is well known 
that bound state formation in < 2 obtains for an arbitrarily weak local potential. 
For higher dimensions, a local potential exceeding a threshold value is required for a 
bound state to form. That such a picture of the bound-state formation applies here is 
supported by simulations on the Hubbard model. In d = (X) |100] a pseudogap is absent, 
whereas a variety of strong-coupling cluster methods all yield a pseudogap [1011 1102] 
without invoking symmetry breaking on a ci = 2 square lattice in the vicinity of half- 
filling. Since L > 2a; is also a signature of a pseudogap (which is mediated by bound- 
state formation), we conclude that dynamical spectral weight transfer also depends on 
the dimensionality of the lattice. The absence of a pseudogap in = oo implies that 
must be some upper critical dimension above which the interactions generated by the 
t/U corrections in Eq. (fT6l) become irrelevant. The precise nature of this fixed point 
remains an open problem. 

A gap in the spectrum is possible only if the single-particle Green function vanishes 
along some surface in momentum space. Along such a surface, the self-energy diverges. 
The imaginary part of the self energy at different temperatures is shown in Fig. (fT3l) . At 
low temperature (T < t^/U), the imaginary part of the self-energy at the non- interacting 
Fermi surface develops a peak at = 0. At T = 0, the peak leads to a divergence and 
hence is consistent with the opening of a pseudogap. As we have pointed out earlier [87]. 
a pseudogap is properly identified by a zero surface (the Luttinger surface) of the single- 
particle Green function. This zero surface is expected to preserve the Luttinger volume 
if the pseudogap lacks particle-hole symmetry as shown in the second of the figures in 

Fig. m- 
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Figure 13. The imaginary part of tlie self energy as the function of temperature for 
n = 0.7. A peak is developed at = at low temperature which is the signature of 
the opening of the pseudogap. The density of states explicitly showing the pseudogap 
is shown in adjacent figure. 



4-2. Mid- Infrared Band 

Naively, doped Mott insulators are expected to either have a far-infrared or an ultra- 
violet or upper-Hubbard-band scale absorption. Hence, one of the true surprises in 
the optical response of the cuprates is the mid-infrared band (MIB). While many 
mechanisms have been proposedjH], no explanation has risen to the fore. A hint as 
to the origin of this band is that the intensity in the MIB increases with doping at 
the expense of spectral weight at high energy and the energy scale for the peak in the 
MIB is the hopping matrix element t. Since the MIB arises from the high-energy scale, 
the current theory which accurately integrates out the high energy degrees of freedom 
should capture this physics. We work in the non-crossing approximation. 



to the Kubo formula for the conductivity where f{uj) is the Fermi distribution function 
and A{u!, k) is the spectral function. In our treatment, the vertex corrections arise solely 
from the interactions with the bosonic degrees of freedom. Since the boson acquires 
dynamics only through electron motion and the leading such term is 0{t^/U'^), the 
treatment here should suffice to provide the leading behaviour of the optical conductivity. 

Shown in Fig. f|T^ is the optical conductivity which peaks at cUmax ~ -St forming 
the MIB. We have subtracted the Drude weight at = to focus sharply on the MIB. 
As the inset indicates, uJmax is an increasing function of the electron filling (n), whereas 
the integrated weight 



Txe^ Jo 

decreases. However, A^efr does not vanish at half-filling indicating that the mechanism 
that causes the mid-IR is evident even in the Mott state. We set the integration 
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Figure 14. Optical conductivity as a function of electron filling, n with the Drude part 
subtracted. The peak in the optical conductivity represents the mid-infrared band. 
Its origin is mobile double occupancy in the lower-Hubbard band. The insets show 
that the energy at which the MIB acquires its maximum value, uj-may, is an increasing 
function of electron filling. Conversely, the integrated weight of the MIB decreases as 
the filling increases. This decrease is compensated with an increased weight at high 
(upper-Hubbard band) energy scale. 

cutoff to f2c = St = 1/m*. The magnitude and filling dependence of f2max are 
all consistent with that of the mid-infrared band in the optical conductivity in the 
cuprates [361 1371 1391 HOl W\\ . We determined what sets the scale for the MIB by studying 
its evolution as a function of U . As is clear from Figure f|T5l) . tUmax is set essentially by 
the hopping matrix element t and depends only weakly on J . The physical processes 
that determine this physics are determined by the coupled boson-Fermi terms in the low- 
energy theory. The Lp\ci^Cii term has a coupling constant of t whereas the Lp\hi scales as 
/U . Together, both terms give rise to a MIB band that scales as cUmaxA — 0.8— 2.21t/f/ 
(see inset of Fig. ([T5|) ). Since t/U ~ for the cuprates, the first term dominates 

and the MIB is determined predominantly by the hopping matrix element t. Within 
the interpretation that Lp represents a bound state between a doubly occupied site and 
a hole, second order perturbation theory with the f^lhi term mediates the process shown 
in Fig. ([5]). It is the resonance between these two states that results in the mid-IR band. 
Interestingly, this resonance persists even at half-filling and hence the non-vanishing of 
A'^efT at half-filling is not evidence that the cuprates are not doped Mott insulators as 
has been recently claimed [8^. 

As the physics in Fig. ([5]) is not present in projective models which prohibit double 
occupancy in the Hubbard basis (not simply the transformed fermion basis of the t- 
J model), it is instructive to see what calculations of the optical conductivity in the 
t — J model reveal. All existing calculations [iOl 1821 1771 183] on the t — J model find 
that the MIB scales as J. In some of these calculations, superconductivity is needed 
to induce an MIB [82] also at an energy scale of J. In others, phonons are invoked 
to overcome the failure of the hard-projected t-J model to yield a mid-infrared band. 
Experimentally [36l l37[ |4T] . it is clear that the MIB is set by the t scale rather than J. 
In fact, since the MIB grows at the expense of spectral weight in the upper-Hubbard 
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Figure 15. Evolution of the optical conductivity for n = 0.9 as U is varied. The inset 
shows the functional form that best describes ojniax- The dominant energy scale is the 
hopping matrix element t since t/U for the cuprates is 0(1/10). 



band, it is not surprising that the t-J model cannot describe this physics as first pointed 
out by Uchida, et aL[37j. The physical mechanism we have identified here, Fig. ([5]) 
clearly derives from the high energy scale, has the correct energy dependence, and hence 
satisfies the key experimental constraints on the origin of the MIB. Since the physics in 
Fig. ([5]) is crucial to the mid-IR, it is not surprising that single-site analysis [83] fail to 
obtain a non-zero intercept in the extreme Mott limit. The non-zero intercept of N^s 
is a consequence of Mottness and appears to be seen experimentally in a wide range of 
cuprates [Ml [M MM HOS] . 

4-3. Dielectric function: Experimental Prediction 

We have shown thus far that there are two branches in the electronic spectral function 
below the chemical potential. Such physics is explained by the dynamical formation of 
a new composite excitation, representing a bound state, consisting of a bound hole and 
a charge 2e boson, (flci^. We demonstrated that for the MIB in the optical conductivity 
such an excitation also appears. Such composite charge excitations should show up in 
response functions which are sensitive to all the charge degrees of freedom, for example, 
the energy loss function, Ql/e{u!,q), where e{u!,q) is the dielectric function. We show 
here that this is the case. 

To this end, we calculate the inverse dielectric function, 

= 7r-J2[ dJ{f{J) - f{uj + ^0) X A{uj + a;', p + q)A(a;', p), 

using the non-crossing approximation discussed earlier. Our results are shown in Fig. f|T6l) 
for n = 0.9 and n = 0.6 for q along the diagonal. Two features are distinct. First, there 
is a broad band (red arrow in Fig. (fT6l) ) with the width of order t that disperses 
with q for both doping levels. It is simply the particle-hole continuum which arises 
from the renormalized bare electron band. The band width is doping dependent as a 



Mottness in High- Temperature Copper- Oxide Superconductors 



35 




Figure 16. The dielectric function, — 31/e(a;,q) for q along the diagonal direction is 
shown for (a) n = 0.9 and (b) n = 0.6. Note only the broad feature indicated by the 
red arrow at n = 0.9 persists aX n — 0.6. 



result of the renormalization of the band with doping. More strikingly, for n = 0.9, 
a sharp peak exists at uj/t ~ .2t. It disperses with q, terminating when q —>■ (vr,7r). 
Physically, the sharp peak represents a quasiparticle excitation of the composite object, 
flcia, the charge 2e boson and a hole. Therefore, we predict that if this new composite 
charge excitation, (pjci^, is a real physical entity, as it seems to be, it will give rise 
to a sharp peak in addition to the particle-hole continuum in the inverse dielectric 
function. Since this function has not been measured at present, our work here represents 
a prediction. Electron-energy loss spectroscopy can be used to measure the inverse 
dielectric function. Our key prediction is that momentum-dependent scattering should 
reveal a sharp peak that appears at low energy in a doped Mott insulator. We have 
checked numerically the weight under the peaks in the inverse dielectric function and 
the sharp peak is important. Hence, the new charge e particle we have identified here 
should be experimentally observable. The two dispersing particle-hole features found 
here are distinct from a similar feature in stripe models [106j . In such models the second 
branch [106] has vanishing weight and whereas in the current theory both features are 
of unit weight. 

4.4- Heat conductivity and heat capacity 

Loram and collaborators [102] have shown from their extensive measurements that the 
heat capacity in the cuprates in the normal state scales as T^. It is a trivial exercise to 
show that such a temperature dependence requires a V-shaped gap density of states as 
a function of energy. The slope of the density of states in the vicinity of the chemical 
potential determines the coefficient of the term. Because the slope of the density of 
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Figure 17. (a) Heat capacity, Cy, and (b) thermal conductivity, r, calculated at 
n — 0.9. The solid lines are a fit to T^. Insert: Density of states for U = lOt are 
evaluated at n = 0.9 and n = 0.7 respectively. 



states decreases as the pseudogap closes, the magnitude of the term should diminish 
as the doping increases. As we showed in the previous section, the boson creates a 
pseudogap. The energy dependence of the gap is shown in the inset of Fig. (|T71) . A 
linear dependence on energy is apparent. We calculated the heat capacity shown in 
Fig.lfTTk) via the relationship C„ = where the internal energy, E, is 



and D{lj) = J^k^i^y^)- As expected, the temperature dependence is quadratic in the 
doping regime where the pseudogap is present as is seen experiment ally [TUT] . As it is 
the boson that underlies the pseudogap, it is the efficient cause of the dependence 
of the heat capacity. In our theory, the steeper slope occurs at smaller doping which 
gives rise to the largest heat capacity at half filling. This doping dependence of the heat 
capacity seems to contradict the experimental observations [107]. A key in determining 
the magnitude of the heat capacity is the spin degrees of freedom. As we have focused 
entirely on the bosonic degree of freedom and not on the contribution from the spin-spin 
interaction terms, we have over-estimated the kinetic energy. Such terms do not affect 
the pseudogap found here though they do change the doping dependence [77]. From Eq. 
( |39|) it is clear that the spin-spin terms renormalize the standard fermionic branch in 
the lower-Hubbard band leaving the new state mediated by ipi untouched. 

Additionally, the thermal conductivity, r(T), can be calculated using the Kubo 
formula in non-crossing approximation. 



The thermal conductivity shown in Fig. (JT711 scales as which is identical to that of the 
heat capacity. However, the system exhibits a larger thermal conductivity as the doping 
increases in contrast to the heat capacity which is decreasing as the doping increases. 
Physically, this signifies that the carriers are more mobile as the doping increases. 




(61) 
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4-5. T -linear Resistivity 

A key theme of this review is that the normal state of doped Mott insulators is dominated 
by dynamical degrees of freedom that could not have been deduced from the UV physics. 
Further, as stated in the introduction, the correct theory of the pseudogap phase should 
also explain the T— linear resistivity. The standard explanation [108] attributes T— 
linear resistivity to quantum criticality. However, one of us has recently shown [T08] 
that under three general assumptions, 1) one-parameter scaling, 2) the critical degrees 
of freedom carry the current and 3) charge is conserved, the resistivity in the quantum 
critical regime takes the universal form. 



As a result, quantum criticality in its present form yields T— linear resistivity (for d=3) 
only if the dynamical exponent satisfies the unphysical constraint 2; < 0. The remedy 
here might be three-fold: 1) some other yet-unknown phenomenon is responsible for 
T— linear resistivity, 2) the charge carriers are non-critical, or 3) the single-parameter 
scaling hypothesis must be relaxed. 

The new dynamical degree of freedom we have identified here fits the bill and 
provides the added ingredient to explain T— linear resistivity. While none of the 
calculations presented here is sufficient to account for the confined dynamics of 
the formation of the pseudogap, the divergence of the electrical resistivity, the v^jcj^ 
feature in the electron operator, and the new feature in the dielectric function all point 
in this direction. Consequently, we assume that v^lcj^ forms a bound state and the 
binding energy is Eb- As a bound state, Eb < 0, where energies are measured relative 
to the chemical potential. Upon increased hole doping, the chemical potential decreases. 
Beyond a critical doping, the chemical potential, crosses the energy of the bound state. 
At the critical value of the doping where Eb = 0, the energy to excite a boson vanishes. 
The critical region is dominated by electron-boson scattering. In metals, it is well- 
known [109] that above the Debye temperature, the resistivity arising from electron- 
phonon scattering is linear in temperature. We make a direct analogy here with the 
electrons scattering off phonons in a metal. Once the boson unbinds, we assume its 
dynamics is purely classical. Since the energy to create a boson vanishes at criticality as 
shown in Fig. (fTSl) . T-linear resistivity obtains. Namely, in the critical region, the energy 
to create a boson vanishes as shown in Fig. (JTSll and hence the resistivity arising from 
electron-boson scattering should be linear in temperature. This mechanism is robust 
(assuming the unbound boson has classical dynamics) as it relies solely on the vanishing 
of the boson energy at criticality and not on the form of the coupling. To the right of 
the quantum critical point, standard electron-electron interactions dominate and Fermi 
liquid behaviour obtains. In this scenario, the quantum critical point coincides with the 
termination of the pseudogap phase, or equivalently with the unbinding of the bosonic 
degrees of freedom. Since it is the bound state of the boson that creates the new charge 
e state giving rise to L/uh > 1 and this state is generated as a result of dynamical 
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Figure 18. Proposed phase diagram for the binding of the holes and bosons that 
result in the formation of the pseudogap phase. Once the binding energy vanishes, 
the energy to excite a boson vanishes. In the critical regime, the dominant scattering 
mechanism is still due to the interaction with the boson. T-lincar resistivity results 
anytime T > uj^, where Wf, is the energy to excite a boson. To the right of the quantum 
critical regime (QCP), the boson is irrelevant and scattering is dominated by electron- 
electron interactions indicative of a Fermi liquid. The QCP signifies the end of the 
binding of fermi and bosonic degrees of freedom that result in the pseudogap phase. 



spectral weight transfer, the T* line defines the temperature below which dynamical 
spectral weight transfer contributes to the low-energy spectral weight. Consequently, 
the mechanism proposed here is experimentally testable. Simply repeat the x-ray K- 
edge experiments presented in Fig. ([3]) below and above the T* line. Above T* the 
integrated weight should be 2x whereas below it should exceed 2x. 

4-6. Towards Superconductivity 

Our emphasis thus far has been on identifying a unifying principle for the normal 
state of the cuprates. As we have seen, strong correlations mediate new composite 
excitations made partly out of the emergent charge 2e boson that results by exactly 
integrating out the high-energy scale in the Hubbard model. An important question 
concerns the relevance of the physics we have identified here to the superconducting 
phase. Equivalently, what role, if any, does dynamical spectral weight transfer play in 
the superconducting state? We answer this question by focusing on a correlate of high- 
temperature superconductivity. As the phase diagram indicates, the superconducting 
region is roughly dome-shaped. Why superconductivity peaks at a particular doping 
level is not known. To offer some insight into this puzzle, we focus on an experimental 
quantity which exhibits an abrupt sign change near optimal doping. As shown in 
Fig. (JT9l) . at a doping level corresponding to the highest superconducting transition 
temperature for a wide range of cuprates, the thermopower vanishes |110j . Consequently, 
the sign change of S occurs at the doping value defining the top of the "dome" . While 
this might be an accident, the fact that the thermopower vanishes at the same doping 
level for most cuprates indicates that the reason might have something to do with the 
superconducting mechanism. 
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Figure 19. a) Universal behavior of the thermoelectric powerpLlO] (290K) as a 
function of planar hole density (Ppi), for various families of hole-doped cuprates. 
All exhibit a sign change at Ppi — 0.23. Above the solid-bold horizontal line, the 
thermpower obeys the functional form, 5'^^°(Ppi) = 392 exp(-19.7Ppi) for 0.01 < 
Ppi < 0.21. Below the solid-bold horizontal line, S'^^°{Ppi) = 40.47 - 163.4Ppi for 
0.21 < Ppi < 0.34. These functional forms were used [110 ] to determine the hole 
doping levels for all the cuprates rather than the widely used empirical formula [lll| 
1 — Tc/T™^^ = 82.6(a; — 0.16)^ which artificially fixes the optimal doping level of 
all cuprates to be 0.16. The thermopower scale is unbiased in this regard and has 
been shown [TTD] to corroroborate independent measures of the doping level even in 
Y123 and Tl-2201 in which it is the oxygen content that determines the doping level, 
b) Maximum transition temperature as a function of the planar hole density using 
the thermopower scale to determine the doping level. Except for three single-layer 
materials, the vanishing of the thermopower coincides with the doping level at which 
the transition temperature is maximized. 



The thermopower measures the thermoelectric voltage induced across a material 
in response to an applied temperature gradient. Microscopically, the thermopower is a 
measure of the entropy per charge carrier. Further, it reveals the nature of the dominant 
charge carriers, being positive for holes and negative for electrons. Should the entropy 
per carrier be identical for particles and holes, the thermopower vanishes. Consider 
hole-doping a Mott insulator. Because transport obtains in the lower-Hubbard band, 
naively a vanishing of the thermopower is expected whenever the number of states above 
and below the chemical potential is equal. In the atomic limit, this corresponds to the 
condition 2x = 1 — x, the solution of which is Xcrit = 1/3- This result is corroborated 
by the large U limit of the thermopower, 

S=-^£,„J^, (63) 

e 1 — a; 
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computed by Beni [112] roughly 20 years before spectral weight transfer was discovered. 
Since 2x and 1 — x are the exact values for the electron addition and removal states, 
respectively, in the atomic limit, it is easy to see that the logarithm is precisely the 
entropy per carrier. The logarithm vanishes at Xcrit = 1/3 which is the exact particle- 
hole symmetric condition for the LHB in the atomic limit. Finite t/U corrections will 
increase L and as a result decrease Xcrit- However, this is not all. The spectral function is 
strongly momentum dependent when t/U ^ 0. As a result, strict particle-hole symmetry 
is not needed to make the thermopower (or even the Hall coefficient) vanish as can be 
seen directly from the exact [113] expression. Both of these effects conspire |1 13] to move 
the doping level at which the thermopower vanishes significantly below the atomic limit 
of X = 1/3. As this change is made entirely from the t/U corrections to the thermopower, 
it is the dynamical spectral weight transfer that is ultimately responsible for the precise 
value of doping at which the thermopower vanishes in the cuprates. Consequently, 
the dynamical spectral weight transfer plays a role in maximizing T^. As a result, 
superconductivity in the cuprates is determined fundamentally by the mixing between 
the high and low energy scales in a doped Mott insulator. The collective degree of 
freedom ip which results from the spectral weight transfer is consequently central to the 
superconducting mechanism. Optical conductivity experiments |114[ I115[ 1116] certainly 
have shown this plainly that the onset of superconductivity results in a decrease in 
the spectral weight in the UHB. It would seem then that the ultimate solution to 
superconductivity hinges on the precise dynamics of the collective charge 2e boson that 
we have shown to exist in the exact low-energy theory of a doped Mott insulator. 
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